Time evolution of damage under variable ranges of load transfer 
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We study the time evolution of damage in a fiber bundle model in which the range of interaction 
of fibers varies through an adjustable stress transfer function recently introduced. We find that the 
lifetime of the material exhibits a crossover from mean field to short range behavior as in the static 
case. Numerical calculations showed that the value at which the transition takes place depends 
on the system's disorder. Finally, we have performed a microscopic analysis of the failure process. 
Our results confirm that the growth dynamics of the largest crack is radically different in the two 
limiting regimes of load transfer during the first stages of breaking. 
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I. INTRODUCTION 

The phenomenon of fracture in heterogeneous mate- 
rials is a complex physical problem which has been of 
great interest to scientists for quite a long time 0, ||, [| . 
A disordered system is understood to be one with ran- 
dom time or/and space dependent breaking properties 
. The random nature of this dependence arises because 
modeling the fracture of heterogeneous materials entails 
dealing with systems made up of many constituents, each 
one having mechanical properties that can be considered 
as being independent of the properties of the other con- 
stituents, but with many body interactions among the 
different parts of the system [Q, g, [| . In heterogeneous 
materials the evolution of the rupture process is radically 
different from the single crack growth mechanism that oc- 
curs in homogeneous materials. Though there is neither a 
complete numerical solution nor analytic solution to the 
fracture problem, we have a better understanding of it 
due to some recent algorithms that have been developed 
to simulate the fracture process j|, |], [| 0, |§] ■ 

Fiber Bundle Models (FBM) form a fundamental class 
of approaches to the fracture problem. They arose in 
close connection with Daniels' and Coleman's seminal 
works on the strength of bundles of textile fibers ^ [l0| , 
and have harbored an intense research activity in recent 
years § 1 1 @ 0, g| g| @ g 0, g, §, f§ 
Fiber Bundle Models are important, despite their very 
simple nature, because they exhibit most of the essen- 
tial aspects of material breakdown. In addition, the 
deep understanding of fracture processes they provide 
has served as a starting point for more complex mod- 
els |n§ g [jj, H H H . In FBMs the material is made 
up of a set of parallel fibers, each having a statistically 
distributed threshold strength or life-time. Besides the 
classical static FBM one can also introduce a dynamic 
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version In the static FBM, the failure process is 

simulated according to a quasi static loading of the mate- 
rial in which the output of the simulation is the ultimate 
strength of the material, i.e, the maximum load above 
which it breaks down. In the dynamic FBM, a constant 
load is maintained on the system and the fibers break by 
fatigue after some time. 

After a fiber breaking, the load acting on it is redis- 
tributed among the intact fibers according to the stress 
redistribution rule associated with the particular load- 
transfer model. There are two standard load transfer 
models comprising the FBM and they correspond to 
the extreme limits of stress redistribution. The global 
load sharing (GLS) redistributes the load of a failed 
fiber equally among the active fibers remaining in the 
system. It is known as the global fiber bundle model 
(GFBM), and assumes long-range interaction among the 
fibers which makes it a mean-field approximation that 
can be solved analytically JU} |l2| . On the other hand, 
the local load sharing (LLSjredistributes the load of a 
failed fiber among the intact fibers that are nearest neigh- 
bors to the failed ones, and is thus known as the local 
fiber bundle model (LFBM). This assumes short-range 
interaction among the fibers and it has not, in general, 
been solved analytically. However, in actual heteroge- 
neous materials the stress redistribution is expected to 
fall in between the two regimes of load transfer. Very re- 
cently, some of us have proposed a load transfer scheme 
with variable range of interaction among the fibers, which 
interpolates between the two extreme cases. Since most 
of the physics of the fracture problem is hidden in the 
stress redistribution, considering within the same model 
the possibility of varying the range of interaction is a 
relatively greater improvement toward a better under- 
standing of the fracture problem. 

Motivated by the results obtained for the static set- 
ting, we have studied a stochastic dynamic fiber bundle 
model, in which fibers break by fatigue over a period of 
time, such that the range of interaction among fibers is 
variable through an adjustable stress transfer function. 
We have observed a crossover from mean-field to short- 
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range behavior for the macroscopic quantities describ- 
ing the fracture dynamics as in the static case. In ad- 
dition, a more detailed inspection of the failure process 
revealed that the microscopic behavior of the system is 
also different as the range of interaction varies. Finally, 
we have also studied the effect of heterogeneity on the 
failure process from macroscopic as well as microscopic 
perspectives. The rest of the paper continues with the 
description of the stochastic model in the next section. 
Section [II is devoted to study the lifetime of the bun- 
dle and the rate at which fibers break when the range of 
interaction changes. The damage spreading for several 
load transfer modalities is addressed in section IV while 



the last section is devoted to conclusions rounding off the 
paper. 

II. THE STOCHASTIC MODEL. 

We assign each fiber to the sites of a square lattice, 
with periodic boundary conditions. Assuming elastic in- 
teraction among the fibers, we may state that the stress 
redistribution obeys a power law [E4| , 
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where Oi nc is the load increment on an intact fiber i at 
a distance r.^ from the failed fiber j, and 7 is a variable 
parameter that controls the effective range of interaction 
among the fibers. 7=0 corresponds to GLS since the 
additional load on each intact fiber as a result of a fiber 
breaking is the same irrespective of its distance from the 
broken fiber. On the other hand, 7 = 00 corresponds 
to LLS since in this limit only the nearest neighbors get 
equal portions of the load of a failed fiber, and <7i„ c =0 
for r > 1. Assuming that there is no dissipation of the 
load of a failed fiber, Eq. (|l|) leads us to the normalized 
stress transfer function, 
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tij being the distance between an active fiber i, with 
coordinates (xi, t/i), and a failed fiber j, with coordinates 
(xj,yj), and A denotes the set of intact fibers. 

One can consider two equivalent approaches to the dy- 
namic FBM (25[ l^j. The time elapsed until the final 
collapse of the system is the lifetime or time to failure of 
the bundle. In the first approach |2^| the lifetime of each 
element is an independent identically distributed random 
variable, i.e, each fiber has a different random lifetime, 
taken from the same statistical distribution (the Weibull 
distribution is a good empirical distribution in materials 
science) and each one breaks if the time elapsed exceeds 
its individual lifetime. This is a quenched model of frac- 
ture where the disorder is introduced once for all at the 
beginning of the process and thus the growth mechanism 
is deterministic. In this version, the effect of the increase 



in stress for a particular fiber i due to the redistribution 
of load from failed fibers is the reduction of its initially 
assigned lifetime U to a new lifetime Ti given by 



U = 



(Tj{t) 
CT 



dt. 



(3) 



where do is the initial stress on fiber i at t = and p 
is the Weibull index, 2 < p < 50. p gives the degree of 
heterogeneity of the system; as p increases the system 
becomes more homogeneous. In this model the next ele- 
ment to break is exactly the one with the lowest lifetime 
at time t. 

, it is considered 
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that in a time step all the intact elements have the same 
mean time to failure. The element that breaks in the 
time interval between two successive failures is selected 
by chance and thus the fiber whose probability of break- 
ing (a function of the load it bears) is the largest is more 
likely to fail. The probabilistic approach is an example 
of the so called annealed disorder since the algorithm is 
stochastic and thus randomness is uncorrelated in time. 
Thus, we start at time t = with N fibers loaded with 
an initial common stress of <Ti(t) — tJi(O) = 1. The mean 
time interval 5 for an individual element to break by fa- 
tigue is, 
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where A is the set of intact fibers and t is the time elapsed 
up to k — 1 breakings, i.e., t = Ep=i &v ^ n the first time 
interval where Ui = 1 Vi, with all the N elements active, 
^1 = jf^p so that the first fiber breaking is completely 

random. This changes with time due to stress transfers 
and fiber failures. The probability of a particular fiber 
breaking j in one sweep of the lattice in the time interval 
S k , is given by 



(5) 



When a fiber breaks, its load is transfered and as a 
consequence, the probability of failure of the receptors is 
increased. In this sense there are no weak or strong fibers, 
all the fibers are equivalent but carrying different loads 
(except for mean field approaches) and thus with differ- 
ent breaking probabilities. It is worth stressing that both 
approaches are equivalent. This can be intuitively un- 
derstood by noting that in both models the load history 
plays a key role. For the quenched setting the fiber that 
breaks is that with the lowest lifetime which in its turn 
depends on the load history. Since the individual times 
to failure are reduced each time a fiber receives load from 
a failed element, the more stressed a fiber is, the more 
likely its lifetime is the lowest. On the other hand, for 
the annealed version the magnitude that is modified by 
the load redistribution is the probability of breaking and 
the load history affects the failure process just as in the 
deterministic version: the more stressed a fiber is, the 
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more likely it breaks. Additionally, we note that this al- 
gorithm is the same as the one used in polymer failure 
pj| and in describing dielectric breakdown [^9| with the 
main distinction that here the broken fibers need not to 
be connected to the single growing cluster as in dielectric 
breakdown. 

Now, the load born by the fiber that has just failed 
is redistributed according to Eq. (|l|) such that in a time 

interval At = t k - t k -i = J2 p =i $p ~ S p =i <^p = $k, the 
load increase on the active fiber i reads as 
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cn{tk) = <7j(ifc_i)5(7,ry) +(Ti(tk-i), 



(6) 



where Uj is the load of the element that has failed in the 
time interval At after k fiber failures. After the redistri- 
bution of load the rupture process continues by applying 
again Eq. (0) that will point to the next fiber to break. 
The lifetime of the material is then given by the sum of 
all the N Ss, T f = J2? =1 Si. 

The complete analytic solution to the dynamic FBM 
as defined above in the GLS (7 = 0) case is feasible. The 
lifetime of the material can be computed as 



T f = 
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that includes also the dependence of the lifetime on the 
system size N giving the mean field result 1/p in the 
thermodynamic limit. 

We should note that there is no avalanche here, unlike 
to the static case 24 , since there is no external driving 



on the system. Once the fibers start breaking by fatigue, 
they continue breaking with time until the final collapse 
of the system, but within a time interval 5k, only one 
fiber may break in a single sweep of the lattice. 

We note that both Eqs. (^) and (jH)) assume a power- 
law breaking rule n p (a) = vq{o j <jq) p . This stems from 
the former assumption that the lifetimes satisfy the 
Weibull distribution. An alternative breaking rule could 
be and exponential hazard rate of the form K e (cr) = 
0exp(?y(cr/CTo)) mainly used for thermally activated pro- 
cesses. On the other hand, k p has the same functional 
form as the Charles power law that describes the subcrit- 
ical crack growth induced by stress corrosion in geological 
materials at constant temperature: v = AKf , where v is 
the crack velocity and Kj is the stress intensity factor for 
mode / fracture. Sometimes, n is known as the stress cor- 
rosion index. Moreover, we emphasize that Eq. (Q) is not 
a real measure of time so that a quantitative comparison 
between the results here shown and those from experi- 
ments is not feasible. Additionally, Eq. (|j) is a simple 
form one can consider for the probability of breaking. 
In principle, one could also include more realistic rules. 
Nevertheless, as we shall see, this simple model is very 
rich and might help understand physical effects present 
in real materials. 
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FIG. 1: Lifetime Tf versus effective range of interaction 7. 
The crossover from mean-field to short-range behavior is ob- 
tained at the same transition point as in the static case 



(P = 2). 



III. LIFETIME OF THE BUNDLE 

We simulate the failure process by large scale numer- 
ical simulations. We first explore the behavior of the 
lifetime as a function of the stress transfer range for a 
fixed level of heterogeneity (p — 2). Then we vary the 
Weibull index such that we get a more homogeneous sys- 
tem with smaller times to failure. The results obtained 
for the lifetime are depicted in Fig. |l| for different values 
of the range of interaction among the fibers and several 
system sizes up to N — 6400 fibers. A crossover from 
mean field behavior to a regime dominated by the short 
range interactions among the fibers is clearly appreci- 
ated. Furthermore, the critical value of the parameter 7 C 
defines a region where the results for global load sharing 
models hold beyond 7 = 0, the true mean-field regime 
for which the load of a failed element is shared equally 
among the surviving elements. For 7 < 7 C the material 
behaves macroscopically as for 7 = 0, that is, the lifetime 
is independent on the system size and does not depend 
on the actual value of the exponent 7. It is not a simple 
numerical task to determine accurately the exact value 
of the transition point due to the stochastic nature of the 
model and the fluctuations of the lifetime of the system. 
In fact, the time to failure of the bundle follows a Gaus- 
sian distribution in what concerns its frequency distribu- 
tion. The width of the distribution depends on the level 
of heterogeneity (controlled by the Weibull index) that 
in turn also influences the lifetimes that become shorter 
as we move to high levels of homogeneity. Within this 
numerical uncertainty we have found that 7 C ~ 2. In- 
terestingly, the same value was found to characterize the 
transition of the bundle's ultimate strength from long to 
short range behavior in the static case of the model [53 . 
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FIG. 2: Comparison of the time to failure obtained for differ- 
ent heterogeneity levels p as the range of interaction varies. 
The inset corresponds to p = 5 while the main figure has been 
drawn considering p — 10. Clearly, the crossover behavior is 
still present although the value of the transition point de- 
pends on p. As the system gets more homogeneous, 7 C shifts 
leftwards. 



The influence of the disorder on the crossover behavior 
can be studied by changing the value of p. Fig. shows 
the time to failure of the material as a function of the 
effective range of interaction for several system sizes. We 
observe that the transition is still present but the range 
where the mean field regime applies is reduced. In par- 
ticular, as the system gets more homogeneous j c shifts 
leftwards to smaller values. Additionally, the true local 
load sharing regime appears to be also slightly shifted 
to the left. When the range of interaction 7 falls below 
a second transition point 7 C 2, the lifetime of the bundle 
becomes again independent on the effective interaction 
among the fibers but it is still size dependent. This later 
behavior can be easily understood by noting that for local 
load sharing schemes the time to failure of the system in 
the thermodynamic limit is zero. We have checked that 
in our model this is actually verified, although the drop 
of the lifetime as the system size is increased is slow. On 
the other hand, wherever the GLS regime arises, the time 
to failure of the system is not size dependent for large N 
(see Eq. ©.) 

Another way to characterize the evolution of the frac- 
ture process is to inspect the rate at which fibers fail. We 
expect two different asymptotic regimes. For long range 
interaction, i.e., below "f c the system should behave in a 
mean field manner. This means that damage is gradually 
accumulated in the material up to a point in which the 
load is too high as to be carried by the remaining fibers. 
Only at this stage of the damage process the rate of fiber 
failures will speed up owing to the small values of the 
very last <5s. On the contrary, in the region where short 
range interaction prevails, the system does not accumu- 



FIG. 3: Normalized number of broken fibers as a function 
of time (also normalized to the lifetime of the bundle) . Two 
groups of curves can be clearly distinguished corresponding 
to the long range and short range regimes. The results have 
been obtained for a system of N — 2500 fibers and p — 2. 



late damage uniformly. In this case there appear regions 
within the material in which stress enhancement takes 
place making the fibers along the crack tips to support 
much more load than other active fibers placed far from 
the clusters of broken fibers. Accordingly, the 8s would 
be modified and there would be more breakings for the 
same time interval. In other words, the breakdown of the 
material occurs suddenly for the very localized regime 
where about 50% of the fibers breaks in a time interval 
of the order of 0.12). In Fig. | we have represented the 
evolution with time of the number of broken fibers Nf for 
different values of 7. Two distinct groups of curves cor- 
responding to the extreme cases can be clearly seen. For 
intermediate values of the effective range of interaction, 
the behavior is more like the case of long range interac- 
tion and may correspond to other load sha ring schemes 
such as the hierarchical fiber bundle model |T§[ |30| . 
Consider that the breaking rate of the bundle is defined 

as 



r(t) 



dN f (t) 
dt ' 



(8) 



with iV>(0) = and N f (T f ) = N. Upon approaching 
the complete failure, the breaking rate scales with the 
lifetime of the material as |31| 



r(t) ~ (Tf - t) 



(9) 



where the exponent £ depends on both the range of in- 
teraction and the Weibull index. However, we can again 
identify two limiting groups of curves for the same value 
of p. Figure ^ shows the rate of fiber breaking for several 
load transfer ranges and a Weibull index p = 2. These 
results confirm the behavior observed for the evolution of 
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FIG. 4: Scaling of the breaking rate r(t) (as defined in the 
text) when approaching the time of material's breakdown. 
The values of 7 are, from right to left, 0, 4, 5, and 20. With 
the exception of 7 = 20, the curves have been shifted to the 
right for the sake of clarity. The least square fit to the data 
gives for the exponent the values £(7 = 0) = 0.5 ± 0.02, 
£(7 = 4) = C(7 = 5) = 0.58±0.02 and Ch = 20) = 0.66±0.02. 
The values of N and p are as of Fig H. The inset shows the 
same quantity for 7 = and p = 10. The scaling exponent is 
in this case £ = 0.86 ± 0.02. 

the number of broken fibers, namely that there is a sharp 
increase in the failure rate when approaching the lifetime 
for the case where short-range interaction dominates the 
damage spreading. The fit to the curves gives ( w 1/2 
for 7 < 7c and £ « 2/3 when 7 is in the range where the 
effective interaction among the fibers can be assumed to 
be very localized. Note that the curves for intermediate 
values of 7 and for the GLS regime have been shifted for 
the sake of clarity. The numerical results are not very 
smooth because of fluctuations but the general trend of 
r(t) confirms the validity of relation (||). As to the de- 
pendence of the above results on the heterogeneity level 
we have observed that the less heterogeneous the mate- 
rial is, the sharper the failure acceleration is in all cases. 
The inset in Fig. [| shows the breaking rate for the case 
of long range interaction and p = 10. The higher value 
of C « 0.86 indicates that much of the fibers break in a 
very small time interval close to the lifetime of the mate- 
rial. As the range of interaction gets more localized, the 
exponent £ approaches unity. 



IV. CRACK GROWTH 

A further characterization of what is going on in the 
fracture process can be carried out by focusing on the 
properties of the clusters of broken fibers. Specifically, 
we have monitored the growth of the cracks inside the 
bundle. Cracks are defined as connected clusters or re- 



gions of broken fibers. Here, we consider a coordination 
number q = 8, that is, each fiber has 8 neighbors. Similar 
results are obtained if we take into account only nearest 
neighbors (q = 4). At the very initial stages, regardless 
of the range of interaction among the fibers, the failure 
of fibers can be assumed to be random, that is, the ini- 
tial cracks are randomly nucleated inside the material. 
This situation changes with time for different load trans- 
fer schemes. By studying the growth of the largest crack 
area C m at each time step, one could distinguish the dif- 
ferent mechanisms leading to the rupture of the material 
as the range of interaction varies. It is worth noting that 
this is just a way that allows to discern between differ- 
ent ranges of interaction and levels of heterogeneity. For 
instance, one can consider instead the linear size of the 
largest crack and study the fractal dimension of the crack 
distribution for different p. 

Figure || shows a microscopic aspect of the material 
breakdown with time for the two limiting cases of load 
redistribution and a system consisting of N = 900 fibers 
(p — 2). Initially, it can be seen that cracks are randomly 
nucleated in the material. At later times the individ- 
ual microcracks coalesce thereby causing a jump in the 
largest crack area. For long range interaction, the nu- 
cleation of cracks continues to be random because all the 
fibers carry the same load and thus they break by chance. 
Therefore the largest broken cluster does not grow lin- 
early with the number of broken fibers. In this case, there 
are isolated cracks inside the material that grow essen- 
tially by coalescence when they meet each other. This is 
the reason why sudden jumps in the area of the largest 
crack are observed in the intermediate stages of the dam- 
age spreading. At the end of the process, the material 
has accumulated many of these cracks giving rise to the 
linear crack growth shown in the inset. 

For localized range of interaction, the mechanism of 
damage spreading is radically different. Again, at small 
times, the cracks are randomly nucleated inside the ma- 
terial. However, as time goes on and more fibers get bro- 
ken, the load is redistributed to the fibers located along 
the crack tips provoking the accumulation of stress in 
these fibers and the appearance of regions where fibers 
bear a huge amount of load. It is thus expected that the 
newly broken fibers add to already existing cracks and 
that a dominant crack appears. From this perspective, 
the largest crack area should scale linearly with the num- 
ber of broken fibers, i.e., C m w Nf. This is indeed the 
case for 7 = 20 as represented in Fig. ||. The straight 
line is a linear fit to a time window in which more than 
150 fibers have been broken. The value 0.90 of the slope 
confirms the above picture. Note that in this case the 
number of coalescence events is smaller than for the GLS 
regime and that after one of such events the linear growth 
of the largest crack is recovered. We also note that, up 
to the intermediate stages of the rupture process, the 
largest crack for a given number of failed fibers is much 
higher in the localized case than the global case indicat- 
ing the formation of a (few) dominant crack(s) in the 
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FIG. 5: Growth of the largest crack area C m with increasing 
number of failed fibers Nf for the two extreme load sharing 
rules here illustrated for 7 = (long range interaction) and 
7 = 20 (short range interaction). The bundle consists of 
N — 900 fibers and p = 2. The straight line is a fit to the 
form Cm = aNf with a — 0.90 ± 0.04. The inset shows 
the evolution of C m up to the macroscopic breakdown of the 
system for the same set of parameters. 



later case. Additionally, at the end of the process there 
are no differences between the two extreme load transfer 
schemes since more than a half of the material is already 
broken and is very unlikely to find a region where iso- 
lated cracks can be formed and grow. Thus, at the final 
stage each additional breaking event occurs at the crack 
tips of existing single dominant cracks. Nevertheless, as 
stated in the previous section the rate of fiber failures is 
quite different in both asymptotic regimes. We shall note 
here that the same behavior as for the mean field regime 
is observed for any value of 7 provided it is below the 
transition point 7 C = 2. 

Figure |^ further substantiates our previous arguments 
by showing how the largest crack area varies as a func- 
tion of the number of broken fibers for several levels of 
heterogeneity. While for 7 = the picture is always 
the same (inset), for 7 = 10, that is, in the localized 
regime, the time taken for cracks to become dominant 
decreases with increasing p. Furthermore, when the sys- 
tem is very homogeneous (p = 10) and local interactions 
prevail, a dominant crack which grows until the mate- 
rial collapses is formed almost instantaneously confirm- 
ing that the mechanism of rupture and crack growth for 
homogeneous materials is radically different from that of 
heterogeneous systems. Nevertheless, for the global load 
case, the change in system's homogeneity does not al- 
ter the way dominant cracks appear and grow (see the 
inset, where no changes, apart form statistical fluctua- 
tions, can be observed) . The reason why this happens is 
given by the way the system gets broken. Equation (|5|) 
tell us that the more stressed a fiber is, the more likely it 



FIG. 6: Growth of the largest crack area C m with increasing 
number of failed fibers Nf for the local load sharing regime 
here illustrated for 7 = 10 and several heterogeneity levels. 
The bundle consists of N = 900 fibers. The inset shows the 
evolution of C m for the long range interaction regime (7 = 0) 
and the same levels of heterogeneity. 



breaks. This always holds except for the global load shar- 
ing case, where the fibers share the same amount of load 
and thus all of them have the same probability to break. 
As the system is more homogeneous (larger p), for local 
load sharing schemes, the probability of breaking for the 
same load is higher so that the appearance of a dominant 
crack is enhanced. Therefore, for long range interaction 
there is no correlated crack growth while for short range 
regimes this is precisely the dominant mechanism since 
the first stages of the damage spreading. 



V. CONCLUSIONS 

We have extended the fiber bundle model with variable 
range of interaction between fibers to the dynamic set- 
ting. As for the static version, two very different regimes 
are identified as the exponent of the stress transfer func- 
tion varies. The lifetime of the material for 7 < 7 C does 
not depend on both the system size and 7 qualifying for a 
mean field behavior. On the contrary, for the short range 
regime, the time to failure of the system systematically 
decreases when increasing the size of the bundle. The 
analysis performed also showed that the exact value of 
the crossover point depends on the level of heterogeneity 
of the system. Besides, we investigated how the material 
approaches its point of total breakdown at both sides 
of the crossover point. The crossover from one regime 
to the other also influences the behavior of the rate at 
which fibers break, explicitly manifested in a power law 
divergence as Tf is approached, but with an exponent 
that depends on the range of interaction. This result is 
relevant from a practical point of view since for the local- 
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FIG. 7: Snapshots of the clusters when nearly 40% of the 
system is broken. The bundle consists of 900 fibers and the 
parameters of the model are shown above the figure. Note 
that varying the heterogeneity level for 7 = does not alter 
the way in which nucleation and crack growth take place. For 
localized regimes, as the system gets more homogeneous, the 
dominant crack appears at early stages and drives the whole 
breakdown of the material. 

ized regime the acceleration of the failure process takes 
place at the very final stages of the rupture process. Al- 
though fibers break by fatigue, one by one, they do break 
in very different time intervals according to the range of 
interaction. In this sense, a global load sharing regime 
is safer, since we get more warnings before the material 
breakdown. On the other hand, the precursory activity 
when the range of interaction gets localized is almost ab- 
sent leading to a sudden breakdown of the bundle in a 
very short time interval. 

The numerical exploration of the damage spreading 
mechanisms under different load transfer regimes fur- 
ther supported the results obtained for the lifetime of 
the system. Regardless of the range of interaction, the 
breaking of fibers is a completely random process at the 
initial stages of the failure process. After some time, 
the mechanism of failure propagation radically changes 



when the exponent 7 varies. In the limiting case of global 
load sharing there is no correlated crack growth in the 
system, whereas for the short range regime the damage 
spreading is driven by a dominant crack, and thus, the 
crack growth is strongly correlated with high stress con- 
centration at the fibers located along the perimeter of the 
dominant cluster of broken fibers. These differences are 
clearly appreciated in Fig. ^, where we represent snap- 
shots of the lattice when nearly 40% of the system is 
broken for several values of 7 and p. For the long range 
interaction limit the material's level of heterogeneity does 
not influence the random nucleation and growth of cracks 
and there are no clearly distinguishable dominant cracks. 
This continues to be so until coalescence drives the fur- 
ther breaking of the material. On the contrary, for lo- 
calized regimes (left column in Fig. |?j), when the system 
gets more homogeneous the dominant crack appears at 
early stages and damage spreading is strongly correlated 
resembling a single crack growth mechanism typical of 
homogeneous materials. When the bundle is heteroge- 
neous, crack growth is still correlated, but in this case 
we can identify more than one large and dominant crack. 
Finally, our results suggest that actually there are only 
two limiting cases relevant to experiments. The one in 
which mean field assumptions apply could be of great 
importance since this will allow the extension of known 
analytic results to ranges of interaction beyond 7 = 0. 
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